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Patch potentials arising from the polycrystalline structure ol material samples may contribute 
significantly to measured signals in Casimir force experiments. Most of these experiments are 
performed in the sphere-plane geometry, yet, up to now all analysis of patch effects has been taken 
into account using the proximity force approximation which, in essence, treats the sphere as a plane. 
In this paper we present the exact solution for the electrostatic patch interaction energy in the sphere- 
plane geometry, and derive exact analytical formulas for the electrostatic patch force and minimizing 
potential. We perform numerical simulations to analyze the distance dependence of the minimizing 
potential as a function of patch size, and quantify the sphere-plane patch force for a particular patch 
layout. Once the patch potentials on both surfaces are measured by dedicated experiments our 
formulas can be used to exactly quantify the sphere-plane patch force in the particular experimental 
situation. 

PACS numbers: 31.30.jh, 12.20.-m, 42.50.Ct, 78.20.Ci 



I. INTRODUCTION 

In distinction to what one is taught in introductory 
physics courses, the surfaces of real metals are not 
equipotentials but are rather described by a locally vary- 
ing surface voltage, known simply as patch potentials. 
Patch potentials exist for several reasons. One is that 
the work function of a crystalline structure depends upon 
which crystallographic plane an electron is extracted 
from [I] [2] . Real metal surfaces are typically composed 
of a network of crystallites with random crystallographic 
orientations, thereby giving rise to a nonuniform poten- 
tial over the metal's surface. In addition, surface con- 
tamination by adsorbates is well-known experimentally 
[3H6] and theoretically [7] to lead to changes in the work 
function. Even for monocrystaline surfaces a spatially 
varying potential has been observed [Sj. Patch poten- 
tials have important implications in various experimen- 
tal disciplines, including gravitational measurements on 
elementary charged particles [5] , tests of the general the- 
ory of relativity PHHT^] . ion trapping [TItTR)] . and the 
physics of Rydberg atoms [17] . In this paper we fo- 
cus on the effects that electrostatic patches can have on 
measurements of the Casimir force [TBH25] . 

Most Casimir force measurements to date have been 
performed in the sphere-plane geometry in order to skirt 
alignment difficulties. The electrostatic interaction is 
used to calibrate the system and to determine the abso- 
lute separation between the sphere and the plane. In the 
idealized case of equipotential surfaces, i.e., no patches 
on the surfaces, the exact analytical expression for the 
sphere-plane electrostatic force (the Coulomb force) is 
well known |24| . valid for arbitrary values of the ra- 
tio D/R, where D is the sphere-plane separation and 
R is the radius of the sphere. In typical experiments 
D I R <C 1 , and the exact expression reduces to its prox- 



imity force approximation (PFA). This approximation 
replaces the sphere by infinitesimal planar surface ele- 
ments and computes the electrostatic force by adding 
plane-plane contributions, as if they were independent. 
For non-equipotential surfaces PFA has also been used to 
compute the electrostatic patch force between the sphere 
and the plane jTSJ HH1 1221 H3]- A further assumption 
in the computation of patch effects has been the er- 
godic hypothesis, that assumes that the actual realiza- 
tion of patches on both surfaces can be well represented 
by statistical properties of their sizes, shapes and volt- 
ages. Since the sphere has a compact cross section the 
sphere-plane interaction can be characterized by an ef- 
fective area of interaction, and the ergodic hypothesis is 
expected to be satisfied when there are many patches 
within the interaction area, thereby providing a fair rep- 
resentation of the patches' statistical properties. 

The main goal of this paper is to derive the exact 
analytical expression for the sphere-plane electrostatic 
patch force which to the best of our knowledge has not 
yet been reported in the literature. Previous theoret- 
ical [HO [22 [22] and experimental [23 [22] works have 
used PFA to address the implications of patch potentials 
on the electrostatic calibration process for Casimir force 
measurements. In this paper we also address the same 
issue using our exact expression for the sphere-plane elec- 
trostatic patch force. 

Measurements of the Casimir force between vacuum- 
separated bodies rely on an electrostatic calibration 
which is performed by sweeping through various values of 
an externally applied potential between the bodies. This 
procedure generates a parabolic force curve as a func- 
tion of applied voltage, the minimum of which identifies 
the minimizing potential. In the absence of patches, the 
minimizing potential is independent of D, and an exter- 
nal voltage equal to the minimizing potential allows for 



the nullification of all electrostatic forces. However, when 
patches are present this is no longer the case. As we will 
show below using our exact expression for the sphere- 
plane patch force, the presence of patches on the samples 
implies that there is no external voltage that can nullify 
the electrostatic interaction - at most the force can be 
minimized by applying a voltage equal to the minimizing 
potential. Also, the existence of patches implies that this 
minimizing potential is, in general, a function of sphere- 
plane distance D. However, we will show that when many 
patches are contained within the effective area of inter- 
action the spatial variation of the minimizing potential 
is suppressed. In this regime the minimizing potential 
may appear, in experiment, to be distance independent. 
Even in this situation the residual electrostatic force is 
non-zero. 

With a detailed knowledge of the patch potentials on 
the sphere and plane surfaces, to be provided by dedi- 
cated ongoing and future measurements, the results con- 
tained here can be used to exactly quantify the con- 
tribution of electrostatic patches to measured signals in 
Casimir force experiments. 

II. EXACT SPHERE-PLANE PATCH FORCE 

We seek the solution to the boundary value problem for 
the electrostatic potential V(x) in the space between the 
sphere and the plane. The potential satisfies the Laplace 
equation 

VV(x) = 0, (1) 

subject to the following boundary conditions: 

V(x)| x£P = V p (x) , ^(x)| xes = V B (x). (2) 

Here V p (x) and V s (x) are the potentials on the plane 
and the sphere, respectively, and P and S denote the 
set of points belonging to the plane and the sphere, re- 
spectively. Once we have the solution for the potential 
we can calculate electric field, E = —W, and thus the 
electrostatic energy in the sphere-plane configuration 

E sp = y^ 3 x(W) 2 , (3) 

where V denotes the volume between the sphere and the 
plane, and e D is the permittivity of vacuum. Finally, 
the electrostatic sphere-plane force is obtained by taking 
minus the gradient of the energy, F sp = — VE sp . 

In the following subsections we will outline the tech- 
niques used to arrive at the exact solution for the poten- 
tial given general electrostatic patchy boundary condi- 
tions, the resulting interaction energy and force, and fi- 
nally the minimizing potential. To this end we will make 
use of bispherical coordinates [57]. The key advantages 
of the bispherical coordinate system for this problem are 
that the Laplacian separates and that the two surfaces on 
which we define our boundary conditions are described 
by 77-constant surfaces. 
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FIG. 1. Slice of the bispherical coordinate system along the z- 
axis. The solid lines are 77-constant surfaces, and the dashed 
lines represent ^-constant surfaces. To span R 3 the coordi- 
nate system pictured above is rotated about the 2-axis; the 
77-constant (77 7^ 0) surfaces will sweep out spheres and the £- 
constant surfaces will sweep out "apples" for < £ < tt/2 and 
"lemons" for n/2 < £ < n. The cross section of the 77 = A 
surface corresponds to the sphere, and 77 = is the plane. 
The relations between the geometrical parameters D and R 
and the bispherical parameters a and A are indicated in the 
figure. 



A. Electrostatic patch energy - Bispherical 
coordinates 

Bispherical coordinates (77, £, (j>) can be used to label 
each point in M 3 . The correspondence with Cartesian 
coordinates is given by the following relations: 

a sin £ cos <f> a sin £ sin (f> a sinh 77 

cosh 77 — cos £ ' ^ cosh 77 — cos £ ' cosh 77 — cos £ 

(4) 

See Fig. [T] for a visual representation. For the sphere- 
plane geometry the adoption of bispherical coordinates 
leads to a significant simplification of the expression of 
the energy through the use of Gauss' divergence theorem 
in curvilinear coordinates 

E sp - ^jf da- Efa = (),£, 0)^(60) (5) 

-^^da-E(77 = A,£,0)I4(£,0), 

where E v is the 77 component of the electric field and 
da = y/H(rf)d^d(f>fj is the oriented measure for the inte- 
gration over an 77-constant surface. The factor yjY,(rj) = 
(a 2 sin £)/ (cosh 77 — cos£) 2 is the square root of the deter- 
minant of the metric induced on an 77-constant surface. 
In these coordinates it is useful to relate A (the 77 coor- 
dinate for the sphere) and a (the location of the foci of 
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the bispherical coordinate system) to the radius of the 
sphere R and the sphere-plane separation D: 

cosh A = l + D/R, 

a = a/OD + R) 2 - R 2 = i? sinh A. (6) 
The Laplacian of the potential in these coordinates is 



v 2 v= ( cosh? ?- cos £)' 



a 2 sin £ 
d ( dV 



d ( dv 
g(v,0 d 2 v 

sin 2 £ 50 2 



(7) 



where g(?y,£) = sin £/ (cosh 77 — cos£). By making the 
ansatz V = VcoshT] — cos £H (r?)H(^)^ ((/>) the Laplace 
equation separates, 



sin 2 ^ ■ 



1 



ff(tj) *(0) 



sin 2 £ 



S'(£) 

+ sm£cos£^y 



(8) 



To arrive at this equation we have divided through by 
csc 2 (£)(cosh?? - cos£)- 5 / 2 #(j7)H(£)*(0). Each of the 
functions 5(f), H(r)), and ^{tf)) can be solved for by sep- 
aration of variables giving the general solution: 



V(t7, £, 4>) — \J cosh rj — cos £ 

oo / 



EE 



Z=0 fc=-i 



(9) 



where A; = I + 1/2, and Afc; and Bki are constants to be 
determined by the boundary conditions on the potential. 
The function is the associated Legendre polynomial 

m- 

The expansion coefficients are obtained by using the 
orthonormality properties of the eigenf unctions which 
serve as the basis for V(rj, £ , (f>). By imposing the bound- 
ary conditions ^ , we find the exact solution for the elec- 
trostatic potential with general patchy boundary condi- 
tions on the sphere and the plane, namely 



V(V,^4>) -EE ^(~ 1 ) fe / W / d£'sin£V fe ^V*(cos£)p- fe (cos£') 



/ cosh 7] - cos £ sinh A^ i] , , / cosh r\ - cos £ sinh A £ (77 - A) , , 
cosh A - cos £' sinhA^A ^ ^ > \\ l-cos£' sinhA^A 



(10) 



Using ( 10 ) the general expression for the electrostatic energy of the sphere-plane system can be written by noting that 
the ^-component of the electric field is E v 



-a 1 (cosh 17 — cos£) J^V(r], £, </>), where we have used the expression for 



the gradient in curvilinear coordinates. By combining the result for E v with the expression for the energy ^ we find 
the general expression for the electrostatic energy 



E sp = yi?sinhA^ ^(-!)* f dn [ ^'e lM ^ ' ) ^(cos£)P £ - fc (cos£') 

— n n f n J J 



= fc=- 

v s (n)v s (n' 



Xf coth Xt A 



sinh A 



v p (n)v p (n') 



x/ cosh A — cos £\/cosh A — cos £' 

v B (n)v p {n') Xe 

-^/coshA — cos£\/l — cos£' sinh A^ A ^/l — cos£^/coshA — cos£' sinhA^A 



2(cosh A — cos£) / ^/l - cosfVl — cos£ 
V P (Q)V S (Q') Xi 



Xt coth A^A 



(11) 



Here J dft is the integration over the "solid angle" 
Jo ^£ Jo* d<f> sin. an d V(fl) = V(£, (f)) in bispherical co- 



ordinates. 



Self-energy terms, both for the sphere and for the 
plane, are contained within Eq.(ll), and their expres- 



sion in bispherical coordinates can be found by taking 
the limit of infinite sphere-plane separation D — > 00. We 
describe the calculation of these self-energy terms in Ap- 
pendix A. In order to obtain the sphere-plane electro- 
static interaction energy these self-energy terms must be 
subtracted from the above expression for the energy. 



B. Electrostatic patch force 

We can find the sphere-plane patch force by taking the 
derivative of the interaction energy with respect to the 
separation, F sp = —(d/dD)E 1 ^ p t . The expression for the 



energy in bispherical coordinates (111 is, however, not 



very transparent for computing the force. This is due to 
the fact that in bispherical coordinates the sphere and 
plane potentials, which prescribe our boundary condi- 
tions, are a function of the sphere-plane separation [29] - 
To circumvent this complication and to connect with the 
natural basis in which the patch potentials are to be mea- 
sured, we transform to the natural coordinate system for 
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each body. Therefore, we make the following change of 
variables (see Appendix A for details of coordinate trans- 
formations from bispherical to spherical or polar coordi- 
nates) 



/ 



sinh A 

dn v B (n)(...) -> I dn s v B (n B ) , * ; —(...), 



(cosh A + cos#) 2 



dn v p (n){. 



dn p v p (n p ) 



Aa 2 



(p 2 + a 2 ) 



r(-), (12) 



where O s = (9, 4>) are spherical coordinates on the sphere 
and Q p = (p, <fi) are polar coordinates on the plane. The 

integration measure J dQ s is given by d(f> d9 sin 9 
with 9 defined as the polar angle on the sphere, and 
J dilp = f Q n d(f> J °° dp p where p is the radius for a po- 
lar coordinate system defined on the plane. We can then 
express the sphere-plane electrostatic interaction energy 
in terms of the natural basis for the two bodies 



E s 



E 



c?fL 



dn b v a (n a )£ a , b (n a ; n b; D)v b (n b ) 



where 
plane 



(13) 

all dependence of the energy on the sphere- 
separation is now contained in the kernels 
£a.b(Qa', n b ; D). Since the functions £ a ,b(n a ', n b ; D) are 
complicated we will place their explicit derivations and 
expressions in Appendix B. By taking the derivative of 
the energy Eq.(13l with respect to D, the exact electro- 



static patch force between the sphere and the plane can 
now be computed 

F sp = [ dn a J dCl b V a (Q a )^ b (il a -n b -,D)V b (n b ), 

a,b— s,p 

(14) 

where F a , b (n a ;n b ; D) = -(d / dD)£ a>b (Sl a] fl b ; D). The 
force above is general for arbitrary boundary conditions 
on the sphere and the plane. 

It is important to emphasize that in these expressions 
the origin of polar coordinate system on the plane is as- 
sumed to be right below the sphere, at the point of closest 
approach between the two bodies. Generally, the mea- 
sured electrostatic patch potential distribution on the 
plane will be done with respect to a different coordinate 
system, say a Cartesian system on the plane. In this 



case, the appropriate change of coordinate system on the 
plane must be performed prior to using Eqs.( |ll|14| ). 

In Appendix C we show how to obtain from our ex- 
act expressions for the energy Eq.(ll) and force Eq.(14) 



for patchy boundary conditions the corresponding well- 
known formulae for the special case of equipotential sur- 
faces. 



Minimizing potential and residual electrostatic 
force 



As mentioned above, in most Casimir force measure- 
ments an external voltage Vo is applied between the two 
surfaces to perform the calibration of the system. By 
sweeping Vo between positive and negative values, the 
total interaction force (or its gradient) versus Vo is mea- 
sured for fixed sphere-plane separation D, resulting in 
force vs. potential "parabola" due to the quadratic de- 
pendence of F sp on Vq. These measurements are then 
repeated for each separation. The minima of each of the 
parabolas defines the minimizing potential, namely 



dF sp 
dV 



(15) 



V =V mi „ 



An explicit exact expression for the minimizing potential 
can be found in this way (see also [THJ [H] for a similar 
approach using PFA). To do so we replace in Eq.(|T4"|) the 
patchy potential on the sphere V s by Vo + V^(SlTJ7 i.e., 
by the addition of the constant applied potential and 
the non-constant patchy one. Alternatively, we can do a 
similar replacement for the patchy potential on the plane 
and, of course, an identical minimizing potential is ob- 
tained. The surface potentials on each of the objects will 
be assumed to fluctuate around the potential given by 
the average work function of the surface material. For 
convenience we will write the potentials in terms of their 
average value and a term describing fluctuations around 
zero. For example, the potential on the sphere becomes 
V s (n s ) = V S + AV s (n s ) where V^ denotes the average po- 
tential and AV^fis) its fluctuations. By solving Eq.(15) 
we find 



V min (D) 



J dn s J dn' s (v s + Av B (n B ))r B , B (n B ; n' s ; d) + j dn p J dn s (v p + Av^rg) Jp, s (n P ; » s ; d) 

fdn s fdn> s F s , s (ci s ;n> s ;D) 
J dn s j dn' s Av a (n s )F s , s (n s ;n' a ;D) + Jdn p j dct s Av p (n p )^ p>s (n p -n s ;D) 



-v s + v„ 



J dn s J dSV B F B , B (Sl B ;Sl B ;D) 



(16) 



In the second line above the integrations against the average potentials have been done. This result can be easily 
understood for the case when there are no patches altogether, AV S — AV p = 0. In this case the minimizing potential 
is the applied potential necessary to nullify the force between an equipotential sphere and plane (see (C9)), and does 
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not depend on D. The integrals against the fluctuating parts of the potential above, AV m in(-D), represent a sort 
of weighted average of the patchy part of the potential and are responsible for all of the spatial dependence of the 
minimizing potential. 

In general, the minimizing potential depends on the sphere-plane separation D through the distance dependency of 
the kernels J- S _ s and J-p, s . This will be shown explicitly in our numerical examples below, where we will also discuss 
the special conditions under which V m i n may appear to be distance-independent even in the presence of patches. 
One should also note that, in general, setting the applied potential Vq equal to the minimizing potential V m - ln does 
not nullify the electrostatic patch contribution to the total sphere-plane force. This can be seen by evaluating the 
electrostatic force at Vq = V m 



* mm- 



sp (v = v min ) = J2 j dQa j dn b v a (n a )F a , b (n a ; n b; D)v b (n b ) - v min J <m s J dn' B F s , s (n s ; n' s ; d) 

a,fe— s,p 

2 J dSlaJ dn b AV a (n a )F a , b (n a ;V b ;D)AV b (n b )- AV^(D) J dn s J dn' s F S!a (n s ;n> a ;D), 



a,fe— s,p * 



(17) 



which is generally different from zero. The second line shows that the electrostatic calibration completely eliminates 
the equipotential component of the force, but does not eliminate the fluctuating part. This residual electrostatic force, 
together with any other voltage-independent interactions (such as the Casimir force), make up the signal in Casimir 
force measurements. 



The expression (17) gives the minimum magnitude that the sphere-plane electrostatic force can take for arbitrary 



surface potentials AV a . Let us now consider what particular form AV a must take in order to minimize the residual 



force given by (17). To do this we will take the variational derivative of F sp (Vb = Vmin) with respect to the surface 



potentials and set the result to zero: 

6F sp (V = Vmin) „ f , . x- . tr > . n. . m at/ m.\ o at/ . m\ f f ac>> t (n . rv . ^AV m{n {D) 



= = 2 Y, I dn b T a Mn a ;Q b ;D)AV b (n b )-2AV min (D) f dtt s f dQ' s j; s (fi s ; tf,; D)- 

dAV a {ila) C^„J J J dAV a (\l a ) 

j m s fdn' s F s , b (n s ]n b ;D)F Sta (M s -,n a -,Dy 



b=s,p 1 



b—s 7 p 



dn, 



F a , b {n a -Q b -D) 



J dn s fdn> s F s , s (n s ;n> s ;D) 



AV b (n b ). (18) 



Since the above equation must be satisfied for all D and the distance dependence cannot be factored out of the 
integral, the function in square brackets can take on nearly any value at each point in the integration domain. This 
leads to the solution for the above integral equation AV a (£l a ) = for which it is easy to verify the F sp (Vb = Vmin) = 0. 
Since this particular choice for the fluctuating potentials gives the extremum for the electrostatic force, any spatial 
variation of the potentials on the surfaces leads to a non-vanishing patch force. This proves that it is impossible to 
nullify the electrostatic force by an externally applied potential when patches are present. 

I 



D. Insights on the spatial dependence of the 
minimizing potential from the PFA 

Before discussing the results of our numerical simula- 
tions, we would like to give some theoretical arguments 
related to the conditions under which the minimizing po- 
tential depends on distance, and to what one can infer 
about patches in the cases where the minimizing poten- 
tial is distance independent. To begin let us explain the 
reasons why an externally applied potential is necessary. 

In Casimir force experiments there exists an intrin- 
sic potential difference between the samples, the contact 
potential, Vc 0n , which is an average surface potential dif- 
ference whose physical origin is the electrical connections 
between the two bodies, differences in work function be- 
tween the samples, and the presence of patches. In the 
following considerations we will neglect the effects from 
connecting wires, solder joints, etc., and focus entirely on 



the contact potential difference arising from work func- 
tion differences and patches. The contact potential leads 
to electrostatic forces between the bodies which can dom- 
inate over the desired Casimir force signal. By applying 
an appropriate bias voltage (the minimizing potential) 
the additional force arising from V con can be nullified. 
Thus, we can understand the nature of the minimizing 
potential via its direct relationship to Vcon. We should 
emphasize that this bias voltage does not nullify the to- 
tal electrostatic force, that has components arising from 
fluctuating patch voltages which are not accounted for in 

Vcon • 

Our goal now is to try to understand the spatial de- 
pendence of the contact potential. To do this we will lay 
out two ideas; the first quantifies the electrostatic sphere- 
plane force via an effective area of interaction, and the 
second relates the contact potential to a weighted aver- 
age of patch voltages. To begin, consider two patchy sur- 
faces interacting within the PFA limit, where the sphere 
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is treated as a large but finite plane. In distinction to 
the case of two infinite planes the sphere-plane force can 
be characterized by an effective area of interaction. This 
is not surprising given that the sphere has a finite cross 
section. Within the PFA we can estimate this effective 
area of interaction by equating the sphere-plane force 
to the product of the plane-plane pressure and the ef- 



fective area of interaction, F s 



Ppp^off, where F sp is 



the sphere-plane force and P pp is the plane-plane electro- 
static pressure. In the PFA limit the sphere-plane force 
can be approximated by the plane-plane energy per area 
E pp , F sp s» 2-7ri?E p p. By noting that E p 
solve for the effective area of interaction giving 



DPpp we can 



A eS ~ 2tt RD. 



(19) 



Now we will connect the idea of an effective area of 
interaction with the contact potential. To do so we will 
assume that each of the bodies are polycrystalline struc- 
tures for which the work function varies above the surface 
due to the different local grain crystallographic orienta- 
tion. Given the variation of the potential over the surface 
we will roughly identify the contact potential with the av- 
erage work function difference between the samples ob- 
served within the effective area. Thus, we can write the 
contact potential formally as 



JV 

^con ^ ^ Vij IVij : 



(20) 



where the indices i and j are used to label the patches on 
the surfaces, N is the number of patches within the effec- 
tive area, and is a normalized weight which is meant 
to roughly account for the fact that patches far from the 
point of the sphere and plane's nearest approach should 
contribute less to the value of the contact potential. 

In order to make these assertions more precise we will 
derive an expression for the minimizing potential (which 
equals the contact potential) within the PFA limit. We 
should stress that this analysis will only roughly char- 
acterize the spatial variation of the minimizing potential 
and is used here to gain physical insight. 

Our starting point is the expression for the sphere- 
plane force found in [T§] : 



sp 



£o 
2 



2ji 



dp p 



(Yo - v(Qp)) 2 

(D + p 2 /2R) 2 



(21) 



where V(f2 p ) is the spatially varying potential differ- 
ence between the two plates, measured at the the po- 
lar coordinate on the plane described by Q p (note that 
V(O p ) should not be confused with the patch potential 
on the plane, Vp(fl p )). By using Eq.(15l we can com- 



pute the minimizing potential (this equation also ap- 
pears in [E5] see Eq.(22)). After replacing V(fl p ) with 
its average value and its fluctuating component, i.e., 



V(ftp) = V + Ay(fip), we find 

r 27T ,, rR , AV(S)p) 

Jo #J d PP [D+^J§Ry 



V 



IT So d P P (D+J/2R)' 



(22) 



In order to proceed we will assume that AV(ft p ) is 
a piecewise constant function, and we will prescribe the 
geometry of the patches. The simplest choice for the 
patch layout is to break each surface into rings centered 
at the point of closest approach and then to divide each 
annulus into several pieces (see Fig. [2j;). Adopting this 
prescription, evaluating the integral in the denominator 
above, and using the fact that D <C R, we can write the 
minimizing potential as 



V 



M K(i) 

EE 

i=l j=l 



<t>i(.i 



dp pDAVij 



2tt Jp.^ R (D + p 2 /2R) 2 ' 
(23) 

Above, the indices i and j label each of the patches, 
i denoting the ring and j denoting the angular sector. 
The term AVij gives the fluctuating part of the poten- 
tial which is constant within the boundary of each patch. 
The total number of rings is given by M, K(i) is the 
number patches in the ith ring, the coordinates <j>j (i) de- 
note the angular boundaries of the patch sectors in the 
ith ring, and pi denotes the outer radius of the ith ring. 
For simplicity we will consider the case where pi = ir avc 
and divide the ith ring into 2i — 1 equal pieces. In this 
way each patch has a fixed area 7rr ave and our expres- 
sion above for the minimizing potential can be simplified. 
This discretization of the surface certainly will not corre- 
spond precisely with a realistic distribution of patches on 
the sample surfaces and is not unique either. Despite this 
limitation we believe that this crude approximation can 
provide some insights into the behavior of the minimizing 
potential. 

After applying this discretization we find the simplified 
expression for the minimizing potential 



F min = V 



M 2i-l 

EE 



2DR L * V V 



(l + (*-l) 2 Si)( 



1 + i 2 



2DR) 

where M = R/r ave (since M must be an integer one 
should take M to be the floor of i?/r ave in numerical 
computations). Since the minimizing and contact poten- 
tials are equal we can use (24) with (20) to obtain an 
expression for the weights 



W.j 



2DR 



i + (*-i) 2 Sb) (i 



' 2DR J 



(25) 



It is interesting to note that these weights arc parameter- 
ized by r ave /2DR, which is the ratio of the patch area to 
the effective area of interaction. One can see that when 
many patches fit inside the effective area of interaction, 
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i.e., rl vc /2DR <C 1, then w l 



oc r" 



J2DR for small % 



(innermost rings), and becomes successively smaller for 

larger rings. At the outermost ring Wij ~ ^Wk^j^- 
Therefore, in this case, all rings contribute to the av- 
erage but their influence is suppressed for large i as l/i 4 . 
In comparison, when the patch area is much larger than 
the effective area of interaction, i.e. r"^ vc /2DR ^> 1, the 
weight for the innermost ring is close to 1, and the weights 
for all other rings are strongly suppressed, roughly pro- 
portional to 2DR/r^ vc . Therefore, in the large patch 
scenario only the patch located at the position of closest 
approach contributes to the contact potential. 

So far we have considered the contact potential for a 
fixed micro-realization of patches, meaning that the volt- 
ages and geometry of each patch have been assigned and 
fixed. This is the case one would encounter in an ex- 
periment. At this stage our calculation cannot proceed 
without a direct knowledge of the patch layout on the 
surfaces. However, based on the simple arguments made 
above it is not hard to make some qualitative statements 
about the expected sample-to-sample fluctuations of the 
contact potential as a function of separation. These ar- 
guments would apply to statistics on an ensemble of min- 
imizing potential measurements performed with samples 
fabricated in the same way. Since the contact potential 
is roughly an average, the more patches which contribute 
to the average the more suppressed will be the sample- 
to-sample fluctuations. This can be roughly understood 
because the uncertainty in an average value scales like 
l/y/N where TV is the number of data points used to 
compute the average. Therefore, a small patch size will 
lead to small minimizing potential fluctuations because 
more patches will determine the contact potential. The 
converse is true for large patches. For the same reasons 
large separations will lead to suppressed minimizing po- 
tential fluctuations since the effective area of interaction 
grows with separation. Likewise, at small separations as 
the ratio r^ ve /2Di? becomes large one should expect large 
sample-to-sample fluctuations of the contact potential. 

We can make these qualitative assertions more con- 
crete by making statistical assumptions about the patch 
voltages. If each of the potentials is assumed to be as- 
signed randomly and statistically independently of one 
another, and with the same variance, then the expected 
variation of the sum on j takes the form Y^/=i AVy — 
±^/2i — 114ms, where Vr ms characterizes the expected 
rms fluctuations of the assigned potential for a single 
patch. This leads to the expected range of minimizing 
potential values given by 



A I 



v ± v ms 



2DR 



£T(i + (*-i) 2 £m)(i + * 2 ^) 



(26) 




FIG. 2. (a) Illustration of the tesselation used for our simu- 
lations; (b) Rings chosen on the sphere; (c) Patch layout on 
the plane. In the numerics the parameter p max is chosen to 
be larger than 5-R (although in (a) p ma x ~ R for illustration 
purposes). For the PFA analysis of the spatial dependence 
of the minimizing potential both plates are tesselated as is 
depicted in (c). 



Thus, in a given experiment, i.e. one micro-realization 
of patches, we would expect the minimizing potential to 



vary in position within a few 0V min (-D) of V. The vari- 
ation is suppressed at large distances because averaging 
is performed over larger and larger effective areas. How- 
ever in the case of large patches, r"^ ve /2DR ^> 1, the i — 1 
term in the sum dominates and the expected variation is 
given roughly by V lms . 

Above we have made some rough arguments in order 
to characterize the expected variation of the contact po- 
tential. We should stress that this expected variation 
will tell us nothing about the minimizing potential in a 
single experiment: equation (26) will not yield a predic- 
tion for the minimizing potential function of dis- 
tance. Rather, the arguments above apply to an ensem- 
ble of measurements of the minimizing potential for dif- 
ferent samples prepared in the same way. The expected 
fluctuations will only tell us the envelope within which 
roughly 66% of minimizing potential measurements of 
various samples will lie. 

The above analysis makes it clear that, in principle, the 
minimizing potential always depends on distance when 
patches are present. However, when the expected vari- 
ation of the contact potential is much smaller than the 
uncertainty in the measurements, the minimizing poten- 
tial will appear in practice to be independent of distance. 
This behavior is expected when the typical patch size is 
much smaller than the effective area of interaction, and 
as we will soon show, can be accompanied by a non- 
vanishing patch force. 
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III. NUMERICAL RESULTS 

In this section we will present our numerical results for 
the sphere-plane minimizing potential ( 16 ) and residual 



electrostatic force (17 1 for various patch sizes. 



To compute the minimizing potential we evaluate 
Eq.(16l in several steps. First, the sphere and the plane 
are tesselated (see Fig. [2]), and then each patch is 
independently and randomly assigned a potential from 
the values (5.15, 5.04, 5.10)V, corresponding to the work 
functions of the three principal crystallographic orienta- 
tions of gold [30]. For each micro-realization, defined 
by a given layout of potentials and patch geometries, 
we compute the minimizing potential using (16). Given 



an ensemble of minimizing potential "measurements" we 
perform statistics to recover the mean value and the 
expected variation. For numerical simplicity we keep 
the geometry of the patch layout fixed in all micro- 
realizations. Another point to acknowledge is that we 
assume that the average potential for the sphere and the 
plane are equal, as if they are made from the same ma- 
terial, and therefore V p = V s . 

The tesselation is performed using a ring-like division 
of the sphere and the plane as is depicted in Fig. [2^,. The 
plane is divided into annuli in the same manner described 
above for the PFA treatment of the spatial dependence 
of the minimizing potential. The outer radius of the ith 
annulus, pi, is given by ir ave and therefore the width of 
each annulus is given by r aV e- The ith ring is divided 
into 2i — 1 equal pieces so that all patches have the same 
area. In order to manage the computation time we set the 
patch potentials to the average potential for radii greater 
than Pmax, which is chosen large enough so that all results 
converge. As we mentioned in our PFA analysis above, 
we acknowledge that this type of tesselation of the plane 
does not correspond with any realistic layout of patches 
on the plane. We believe, however, that this approxima- 
tion is sufficient to gain some insight about the patch size 
dependence of the force. Next, the sphere is tesselated in 
a similar manner by being divided into a fixed number 
of latitudinal lines any neighboring pair subtending polar 
angle AO, afterward the lines of longitude dividing a given 
ring are chosen so that all of the patches have the same 



area. We choose AO to take the value RAO 



For 



this particular tesselation the polar angle describing the 
upper boundary of the ith ring is given by 0i = ir — 7 -^ L i 
where we have chosen i = to correspond with the polar 
angle of the "south pole" . With the particular choice we 
have adopted for the polar angles of the rings we find the 
number of patches in the ith ring on the sphere is given 
by |Floor[(2i? 2 /r2 vc )(cos0 l _ 1 - cos0;)]|. 

The numerical results for the minimizing potential are 
presented in Fig. [3j All of the plots employ the same 
convention: the solid line denotes the average minimiz- 
ing potential as a function of distance, and the dashed 
lines are computed from the standard deviation of the 
all minimizing potential values at each separation. Thus, 
by definition, the envelope created by the two dashed 



lines contains roughly 66% of all minimizing potential 
measurements. The data points indicated by the various 
plot markers are the values of the minimizing potential 
from five random micro-realizations. The figure illus- 
trates the spatial dependence of the minimizing poten- 
tial on patch size. In line with our PFA analysis, as the 
patch size grows (from (a) to (c)) the expected minimiz- 
ing potential fluctuations are enhanced. Additionally, for 
each patch size the effective area of interaction decreases 
as the sphere-plane distance becomes smaller. Thus, the 
shorter the distance, the fewer the patches which con- 
tribute to the average that determines the minimizing 
potential, and hence the larger are the fluctuations. 

To obtain the sphere-plane electrostatic patch force 
Eq. ( 14 ) , we employ the tesselation described above, see 
Fig. |2J Once we have an ensemble of realizations we 
find the average patch force and its expected fluctuations. 
The results of our simulations are depicted in Fig. |4j We 
also compare these numerical results with the patch force 
computed by performing an average over voltage micro- 
realizations analytically while keeping the patch geomet- 
rical layout fixed. To do this we compute the two-point 
voltage correlation functions on the plates, where we as- 
sume that the voltages on each patch are statistically in- 
dependent, random variables (so-called "quasi- local cor- 
relation", see Eqs. (10-12) of [53]), 

(V a (x)V b (yf)) v = 5 ab V r 2 ms e i (x)e i (x / ), (27) 

i 

where the function 0i(x) is 1 for points x within the ith 
patch and otherwise. Using this equation in |l4| ) we 
compute the ensemble averaged patch force (F sp ) v . In 
Fig. ([4]) we compare the exact numerics with the patch 
force computed via the voltage correlation method. As 
seen in the figure, the agreement between the two is ex- 
cellent. In the inset wc compare the residual electrostatic 
force ( 17 ) with the above patch force ( 14 ) . The two forces 



are very similar, illustrating that the electrostatic calibra- 
tion does not nullify the patch force. 



IV. CONCLUSIONS 

In this paper we have derived the exact expression for 
the sphere-plane electrostatic potential, interaction en- 
ergy and force for patchy boundary conditions. With 
knowledge of the potentials on the sample surfaces used 
in Casimir force experiments in the sphere-plane geom- 
etry, to be measured by dedicated experiments, these 
expressions can be used to exactly quantify the sphere- 
plane force due to patch potentials. 

As an added benefit of the exact solution we have been 
able to derive an exact formula for the minimizing poten- 
tial used to calibrate many Casimir force experiments. 
We have found that when patches are present the min- 
imizing potential always depends on position. However, 
the magnitude of the spatial fluctuations are controlled 
by the ratio of the typical patch size to the effective area 
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FIG. 3. Minimizing potential as a function of distance for 
various patch sizes: (a) 100 nm, (b) 300 nm, and (c) 600 nm. 
The solid curve is the average minimizing potential computed 
from 400 micro-realizations. The dashed curves enclose the 
range of expected fluctuations computed from the standard 
deviation of all minimizing potential values at a fixed dis- 
tance. The data points denote the minimizing potential for 
five random micro-realizations. The radius of the sphere is 
R = 150/im. 
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FIG. 4. Sphere-plane electrostatic patch force (14 1 as a func- 
tion of separation for different patch sizes: 100 nm (dotted), 
300 nm (dashed), and 600 nm (solid). The data points are 
the average patch force, the error bars are not visible at scale 
shown, and the lines are the patch force computed with the 
voltage correlatio n m ethod (14ms = 45 mV). Inset: Ratio of 
the residual force ( 17 1 to the patch force ( 14 1 as a function of 



separation for different patch sizes, 
is R = 150/^m. 



The radius of the sphere 



patches contribute to the average (by having comparable 
weights) the expected fluctuations of the minimizing po- 
tential are suppressed. This happens when many patches 
fit within the effective area of interaction. In contrast, 
the expected fluctuations are greatly enhanced when the 
typical patch area is comparable to the effective area of 
interaction, as the patch at the point of nearest approach 
dominates the weighted average. 

Furthermore, we have derived an explicit formula 
for the residual electrostatic force (17) that persists in 
Casimir force measurements after an electrostatic cali- 
bration has been performed. By analyzing the residual 
electrostatic force we proved that the electrostatic inter- 
action cannot be nullified by an applied field unless there 
are no patches. This behavior was also verified in our 
simulations. 
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Appendix A: Electrostatic self-energies 

In this Appendix we compute the large separation limit 
of the electrostatic sphere-plane energy (111 in order to 



identify the self-energies of the sphere and the plane in 
bispherical coordinates and show that they correspond 
to the well-known self-energies of each of these objects. 
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We start with the computation of the self-energy for 
the sphere. First, we briefly recall the derivation of its 
electrostatic self-energy in spherical coordinates. The 
electrostatic potential for an isolated sphere of radius R 
with surface patch potentials V s (fl) = V s (9, </>) is given by 
the known solution of the Laplace equation in spherical 
coordinates, 



Xt 



f) 



e+i 



e=o m=- 

dQ' s e im ^-^PP(cosd)Pf m {cose')V B (fl' s ).(Al) 



Similar to the procedure used in the main text, we use 
the Gauss divergence theorem to express the sphere self- 
energy as 



self 



(A2) 



r=R 



where the integration is over the sphere's surface. In this 
way we obtain the known electrostatic self-energy of the 
sphere in spherical coordinates 



dn„ / <m' 



Er lf = |^E E w + 1/2) (-d- 

e=o m =-e 

xe im ^-^PJ n '(cos9)Pi m {cose')V s (n s )V s {n' s ), (A3) 

where we used that Xi — £ + 1/2. 

Next we show that by taking the large distance limit 
in (11), the piece which is a quadratic function of the 



sphere's potential gives exactly this result for the sphere's 
self-energy. When D — > oo, or equivalently when A — > oo, 
we get from (11) 



Ef f = lim Y A ^ + 1 /2) ( _ 1)fc 



e=o k- 



2tt 



dn / do! 



x e ife ^-^)pfc(cos^)P7 fe (cos^)K(^)K(fi') 



(A4) 



We use the following identities that relate bispherical 
coordinates with spherical coordinates centered on the 
sphere 



cos£ 
sin£ d£ 



1 



(cosh A — cos£) n 



1 + cosh A cos 9 
cosh A + cos 9 
sinh A sin 9 
cosh A + cos 9 ' 

sinh 2 A 
(cosh A + cos6>) 2 
(coshA + cosfl)™ 
sinh 2 " A 



(A5) 
(A6) 

sin 6» dO, (A7) 
(A8) 



For A — > oo we obtain cos£ — > cos 9, d£l — > dQ, s , and 
therefore the sphere self-energy in bispherical coordinates 
( A4 ) is identical to that in spherical coordinates ( A3 ) , as 
expected. 



We now consider the self-energy for the plane. The 
electrostatic potential for an isolated plane located at 
z = having surface patch potentials V p (Q) — V p (r) 
(r = (x,y)) is 



V(x) 



d 2 k 



(A9) 



where Vp(k) is the 2D Fourier transform of the patchy 
potential. The self-energy is 

= Y j d2r j rf2l -%( r )/( r - r ')^( r '),(Aio) 

where J(r-r') = (2tt)- 2 / d 2 k|k|e ik -( r - r '' = -(27r)- x |r- 
r'|" 3 - 

As for the case of the self-energy for the sphere, we take 
the D — > oo limit (or, equivalently, the a — > oo limit) of 
(11) and focus on the term that depends on the plane. 



We get 



00 £ \ 2 

p 2 ^ ^ 2tt 

1=0 k=-i 



E E £(-*)* / < m i < m ' 

£=0 k=-i 
y/l — COS^Vl ~ C0s£' 



(All) 



In order to show that this is identical to ( A10 1 we first 
perform the summation over k in Eq.( |All[ ). For this 
we use the addition theorem of Legendre polynomials, 

EL-^(- 1 ) feeife( ^ <A ' ) - P /( cos 0-fT fe (cosf) = P(cos 7 ), 
where 

cos 7 = cos £ cos £' + sin£sin£'cos((/i> — </>'). (A12) 

In the resulting equation we then change coordinate sys- 
tem from bispherical coordinates to polar coordinates on 
the plane (77 = 0), centered at the point of closest ap- 
proach between the sphere and the plane, i.e. right below 
the sphere. The polar coordinate p relates to the bispher- 
ical coordinate £ through p = acot£/2. The following 
identities are useful to relate both coordinate systems: 



cos£ = 
sin£ = 



2(p/q) 
(p/a) 2 + V 

2V2pdp 



sin £d£ 

VI- cos £ = ~ a 2 [1 + (p/a) 2 } 3 / 2 ' 



(A13) 
(A14) 
(A15) 



Then Eq.(All) takes the form 



E s p clf = lim ^| [dn p [ dii'v p (n p )v p {n' p ) 

OO 

x^A 2 P £ (cos 7 ), (A16) 
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where ft p — (p, cj>) are the polar coordinates on the plane, 
and dflp is the corresponding measure. The above sum- 
mation can be evaluated using the generating function of 
Legendre polynomials, that verifies 



S(t,u) = J2t e P e (u) 



(A17) 



V 7 ! - 2tu + t 2 ' 
Recalling that Xg = I + 1/2, we obtain 

± A,Vft(«) =(«!<! + «» + i)5( t .«) 
e=o v 7 

_ 1 - 10t 2 + t 4 +4t(l + i 2 )w 
4(l + t 2 -2iu) 5 / 2 

= 0(t,u). (A18) 

Setting t = 1 we derive the identity ^|-^( cos 7) = 



[2(1 - COS7)]- 3 / 2 . Using (|A13[) and (|A14[) in the def- 



inition of cos 7 and taking the a — > 00 limit we obtain 
1 — cos 7 s» 2|r — r'|/a 2 . Hence ( A16 ) is equal to 



^J d2r J d 2 r'V p (r)\r-r'\- 3 V p (r'), (A19) 



which is exactly the plane self-energy in polar coordinates 
(lAlOj), as it should. 



Appendix B: Derivation and expression of 
sphere-plane energy kernels 

In this Appendix we give the derivation of the kernels 
appearing in Eq.( 13 ). We first write the energy ( 11 ) in a 



compact form, E sp = Eg 



-E s ' p - 



EPf + EPf , where the 



superscripts denote which surface potentials contribute 
to each term. For example, let us consider the E^-term 



El-; = ^R S mhAj2 J2 ^(~ 1 ) fe / dn J dVe^-^ (cos £)Pi k (cos O 



e=o k=-t 

v s (n)v s (n>) 



\J cosh A — cos cosh A — cos £' 



Xe coth Xi A ■ 



sinh L 



2 (cosh A — cos£) 



(Bl) 



In order to simplify the above expression we first notice that the second term in the parenthesis is independent of 
i, which allows us to use the completeness relation for the associated Legendre polynomials 



00 e 



E E ^(~l) k e lk{ ^' ) P e k (cosOPi k (cose) = S(<P-^)6(cos^-cosa, 



=0 k- 



2vr 



(B2) 



and perform the f2 integral. We also notice that the first term in the parenthesis can be written as coth Xi A 



°° '-2X/Ar 



where the prime in the summation symbol indicates that the n = term has to be taken with half 



weight. By using the completeness relation Eq.(B2), the addition theorem for Legendre polynomials, and the identity 
Eq.(C2) we can express Ejj£ as 



sinh 2 A V s 2 (n) e Q 
(cosh A — cos £) 2 2ir 



E 



,— Ar, 



dn / dn 



, aVsi^Vsig) g(e- 2A ",cos7) 
\f cosh A — cos £y/ cosh A — cos £ ' ' 



(B3) 



where we have introduced the function Q(t,u), defined in Eq.(C2). Using the identities in the previous Appendix, 
that relate bispherical and spherical coordinates, we obtain 



r / dn s v s 2 (n s 



2n 



E 



-At. 



, () , , Rsinh 3 AV s (n s )V s (n' s ) g(e- 2A ",cos7 5 ' s ) , . 

' ' s ' (coshA + cos0) 3/2 (coshA + cos#') 3/2 ' 



where cos7 s,s is given by (A12) for the two bispherical coordinates (£,</>) and (£',(/>') written in terms of spherical 
coordinates, namely 



cos 7 



(1 + cosh Acos0)(l + coshAcosfl') + sinh 2 A sin 9 sin 9' cos(0 - </>') 
(cosh A + cos 9) (cosh A + cos 6') 



(B5) 



The same kind of calculations can be performed for the other terms EjjjP, E|£ 8 , and EP^ P in the energy (11 ). In this 
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way we derive the energy kernels that appear in (13): 
S SiS (n s ;{l' s ;D) 



^—SU -^Wcos 6 -cos 0') + — # lim g(l,cos 7 s ' s ) 

4 4-7T A-voo 



self —energy 

iJsinh 3 A e- A,i 0(e- 2Ar \ 



cos 7 



2ir ^ (cosh A + cos 9) 3 / 2 (cosh A + cos 9' ) 



3/2 - 



£ (n-n-m-s fn-n ^)--^V fl 2 sinh 2 a g(e - 2 ^ COS7 s, P) 

S ' Pt S ' P ' j ^ Zs '^' ^ " 7T ^ ( C osh A + COS0) 3 / 2 (i? 2 sinh 2 A + p 2 )3 /2 > 



£0 



4e 



47T X - 



,/|3 



E 



i? 3 sinh J A e 



-An 



-2 An 



COS 7 P ' P ) 



} (i? 2 sinh 2 A + p 2 f/ 2 (R 2 sinh 2 A + p' 2 ) 3 / 2 



(B6) 
(B7) 
(B8) 



self —energy 



The dependency of these kernels on the sphere-plane separation D is encoded in A, which we recall depends on D as 
cosh A = 1 + D/R. The functions cos 7 s,p and cos 7 P:P are obtained by expressing (A12) in terms of the corresponding 
spherical or polar coordinates: 



s,p 



cos 7 



cos7 p,p = 



P 



R 2 sinh A 1 + cosh A cos 9 2Rp sinh A sinh A sin 9 



p 2 + R 2 sinh A cosh A + cos 9 p 2 



R 2 sinh A cosh A + cos 9 



(p 2 - R 2 sinh 2 A){p' 2 - R 2 sinh 2 A) + AR 2 pp' sinh 2 Acos(</> - 4> r ) 
{p 2 + R 2 sinh 2 A)(p' 2 + R 2 sinh 2 A) 



(B9) 
(BIO) 



The first two terms in ( B6 1 and the first term in 



correspond to the self-energies of the sphere and the plane, 
respectively, already derived in the previous Appendix. There should be removed from the energy kernels when 
computing the interaction energy or force. 

I 

Appendix C: Equipotential case 



In this Appendix we show how to obtain from our gen- 
eral formula (11 1 the electrostatic energy when both the 



sphere and the plane are equipotentials. First we will 
consider the component of the energy depending quadrat- 
ically on the sphere's potential 



f^E E ^(-!) fe [ dn I ^ 4fc( ^V(cos£)P^cosO 

n n j /> J J 



e=o k=-t 
a 



-^/cosh A — cos £ 



V^cosh A — cos £' 



Xf coth A* A 



sinh A 



2(cosh A — cos£) 



(CI) 



The ^-integrals can be done directly giving Kronecker identities 
deltas which collapse the sum on k, and the integrals 
over £ can be done after using the identity —^==== = 

\/2e _A / 2 X^n'Lo -Pn(x)e~ An , and the orthonormality prop- 
erties of the Legendre polynomials giving the following 



<%sinffl(cosfl 1 = ^ e ~^ A , 

„ v cosh A — COS £ \g 

() £ ' € e £ (cosh A — cos£) 3 / 2 sinh A 
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After performing these integrations Eq.(Cl| reduces to 



E" q s _ p =2ne aV s 2 e~ 2Af A (coth A<A + 1) 



e=o 



--2ne aV 2 V ^ 

L — J si 



-A.A 



sinh \eh. 



1=0 

Note that the previous equation can be written as 



E s e '*_ p =2ne aV s 2 £ 



sinh^f sinh ^ A 



(C2) 



(C3) 



which after some rearrangement reduces to the well- 
known form for the energy for the prescribed case of 
cquipotcntials [24 



oo 



sinh An 



(C4) 



With the same identities one can derive the other terms 
contributing to the electrostatic energy. The cross terms 
are given by 



E s ' p = E p ' s 

cq— p cq— p 



2ire aV s V p ^ 



e=o 



sinh A<A ' 



(C5) 



and the component depending on the square of the 
plane's potential can be written as 



EP q p „ p =2ne aV 2 coth A ^ A - 



(C6) 



£=0 



The contribution to the electrostatic energy from the 
preceding equation contains a divergence which results 
from the infinite self-energy of the plane. Note from 
Eq.(All| and using Eq.(C2) that for the case of an 



equipotential on the plane the plane self-energy reduces 
to E P ^ pself _ enorgy = 2TTe o aV 2 J2T=0 1 - % subtracting 
this divergence the plane-plane contribution to the elec- 
trostatic energy becomes 



e p,p „ E P:P 

cq— p cq—p, self— energy 



=27re aV 2 J2 e 2X e A_ 1 



--2ns aV 2 Y; 



-\tA 



1=0 



sinh \o A 



(C7) 



By combining all of the terms which contribute to the 
energy after all divergences have been removed we find 
the final form for the equipotential energy 



oo ^ 

E cq _ p = 2ne (V s ~V p ) 2 aJ2—r 



sinh An 



(C8) 



The expression for the force in the equipotential case 
can be derived by taking a derivative with respect to 
sphere-plane separation giving 



_ /Tr T r . , v"-r coth A — ncoth An . 
= ~^o(V s - V p ) 2 £ ■ (09) 
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